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channel  state  infonnation  (CSI)  mismatch,  and  imperfect  power  control.  Numerical  results  show 
that  the  proposed  adaptive  transmission  method  has  much  greater  system  capacity  than 
conventional  non-adaptive  MC  direct-sequence  (DS)-CDMA  system. 

Index  Terms  —  Code  division  multiple  access,  Multiuser  channels,  Fading  channels,  Frequency 
hop  communication,  Resource  management 


i 


I.  INTRODUCTION 


Multicarrier  code-division  multiple  access  (MC-CDMA)  methods  have  been  extensively 
investigated  and  proposed  for  future  mobile  radio  systems  due  to  their  high  spectral  efficiency, 
the  interference  suppression  capability,  the  exploitation  of  frequency  diversity,  and  the  ability  to 
adapt  to  time-varying  channel  conditions  [1,2, 3, 4],  During  the  past  decade,  several  MC-CDMA 
systems  with  adaptive  frequency  hopping  (AFH)  were  proposed  for  cellular  communication 
systems  in  [5,6,7]  and  have  attracted  significant  attention  due  to  their  excellent  spectral 
efficiency.  In  these  systems,  the  data  of  each  user  is  multiplexed  over  one  or  several  substreams, 
and  multiple  subcarriers  are  employed  to  transmit  the  substreams  of  all  users.  Direct  sequence 
spread  spectrum  (DS/SS)  codes  are  assigned  to  all  substreams  in  the  system.  Each  substream  is 
transmitted  on  one  subcarrier.  The  subcarriers  are  allocated  using  the  knowledge  of  the  channel 
state  information  (CSI).  MC-CDMA  with  AFH  exploits  both  frequency  and  multiuser  diversity 
and  was  shown  to  improve  on  non-adaptive  MC  direct-sequence  (DS)-CDMA  system  in  [8]. 

In  [6],  a  MC-CDMA  system  with  AFH  was  proposed  for  the  forward  link.  In  this  system,  one 
substream/user  is  employed,  and  each  user  selects  the  subcarrier  with  the  largest  fading 
amplitude  for  transmission.  This  allocation  algorithm  provides  significant  perfonnance 
improvement  over  other  diversity  techniques  such  as  non-adaptive  MC  DS-CDMA  with  the 
maximum  ratio  combining  (MRC)  [8],  However,  this  simple  allocation  algorithm  does  not  take 
into  account  the  multiple-access  interference  (MAI)  present  in  realistic  wireless  systems.  A  near- 
optimal  allocation  algorithm  was  utilized  in  [9]  to  maximize  the  total  signal  to  interference  and 
noise  ratio  (SINR)  and  was  shown  to  improve  on  the  performance  of  the  simple  allocation  policy 
in  [6]. 

For  the  reverse  link,  a  MC-CDMA  system  with  AFH  was  investigated  in  [5].  In  this  system, 
multiple  substreams  are  employed  for  each  user,  and  random  signature  sequences  are  assigned 
for  all  substreams,  resulting  in  enhanced  MAI  and  intra-user  interference.  A  sub-optimal  water¬ 
filling  (WF)  allocation  algorithm  was  proposed,  and  it  was  demonstrated  this  system  has  better 
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performance  than  DS-CDMA  system  with  RAKE  receiver  and  the  non-adaptive  MC  DS-CDMA 
system  [5].  However,  the  WF  algorithm  offers  limited  protection  to  weaker  users.  Thus,  it  suffers 
from  the  near-far  problem  caused  by  short-term  fading  and  imperfect  power  control.  As  a  result, 
the  bit  error  rate  (BER)  of  the  WF  method  is  degraded  in  moderate-to-high  signal-to-noise  ratio 
(SNR)  region. 

In  this  paper  and  [10,  17],  we  investigate  an  improved  MC-CDMA  system  with  AFH  for  the 
reverse  link.  This  system  exploits  realistic  quasi-synchronous  reverse  link  model  [11]  and 
utilizes  orthogonal  spreading  sequences  for  all  substreams  of  the  same  user,  thus  eliminating 
intra-user  interference.  Moreover,  we  propose  a  new  allocation  algorithm  to  overcome  the 
limitations  of  the  WF  algorithm  in  [5],  and  demonstrate  the  resulting  BER  improvement. 

MAI  severely  limits  perfonnance  of  spread  spectrum  systems.  While  it  has  been  demonstrated 
that  multiuser  detector  (MUD)  can  greatly  improve  spectral  efficiency  of  DS-CDMA  [12]  and 
non-adaptive  MC-CDMA  [1,13],  to  the  best  of  our  knowledge,  it  has  not  been  utilized  in  MC- 
CDMA  with  AFH.  In  this  paper,  we  employ  the  MUDs  jointly  with  the  proposed  allocation 
algorithm  by  utilizing  the  SINR  analysis  [10,33].  Performance  gains  achieved  using  the  MUD 
combined  with  the  proposed  allocation  method  are  demonstrated.  To  reduce  the  complexity  of 
the  allocation  algorithm  combined  with  linear  MUDs,  recursive  update  of  the  matrix  inverse  is 
utilized. 

In  previous  research  on  MC-CDMA  with  AFH  [5,  6],  perfect  knowledge  of  the  CSI  at  the 
transmitter  and  independent  fading  for  all  subcarriers  were  assumed  in  the  performance 
evaluation.  These  assumptions  are  too  optimistic  for  realistic  wireless  channels.  The  impact  of 
noisy  CSI  measurements  was  considered  in  [6,  7].  However,  imperfect  CSI  is  mainly  caused  by 
the  feedback  delay  associated  with  sending  channel  allocation  instructions  from  the  base  station 
to  mobile  users  in  a  rapidly  varying  fading  channel.  While  the  impact  of  imperfect  CSI  has  been 
studied  extensively  for  other  adaptive  transmission  systems  [14,  15,  22,  23],  this  investigation 
has  not  been  carried  out  for  the  MC-CDMA  system  with  AFH.  In  addition,  it  is  well  known  that 
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correlated  subcarriers  reduce  the  frequency  diversity  of  MC  systems  [16],  and  the  assumption  of 
independent  subcarriers  tends  to  overestimate  the  system  performance  [17].  In  this  paper,  both 
imperfect  CSI  and  correlated  subcarriers  are  considered  in  the  evaluation  of  system  perfonnance. 

The  remainder  of  this  paper  is  organized  as  follows.  In  section  II,  we  describe  the  reverse  link 
model  of  the  proposed  MC-CDMA  system  with  AFH  and  MUDs  employed  in  this  paper.  Section 
III  analyzes  the  limitations  of  the  WF  allocation  method,  and  describes  the  proposed  allocation 
algorithm.  Numerical  results  and  conclusion  are  contained  in  sections  IV  and  V,  respectively. 
Symbols  frequently  used  in  this  paper  are  summarized  in  Table  I. 

II.  SYSTEM  MODEL 

A.  Transmitted  Signal 

System  diagram  of  MC-CDMA  with  AFH  is  illustrated  in  Fig.  1 .  In  the  system  with  K  users, 
the  total  bandwidth  W  is  divided  into  M  subcarriers  with  equal  bandwidths  A f=  W/M .  The  data 
stream  of  each  user  is  multiplexed  over  N  substreams,  and  all  substreams  are  spread  by  spreading 
codes  in  time  domain.  The  low-pass  equivalent  signal  at  time  t  for  the  nth  substream  of  the  Mi 
user  is  (see  Table  I) 

00  PG- 1 

xk,n(t)=-\flEh  Z  Z  bk,n(l)Ck,n(i)h(t-lTb-iTc )  (1) 

/=— oo  z— 0 

where  bk,n(l)e  {— 1  ,+l }  are  the  binary  phase  shift  keying  (BPSK)  modulated  infonnation  symbols. 
The  spreading  code  vector  for  this  substream  is  defined  as  Ca,;i=[ca,„(0)  CA,;i(l)  ...  ckn(PG- 1 )] 7 

with  ||Cfc„||=l,  where  ||*||  is  the  Euclidean  nonn,  and  PG  is  the  processing  gain.  The  chip 

00 

waveform  h(t)  satisfies  fh2(t)dt=\. 

-oo 

After  spreading,  the  substream  ( k ,  n)  is  assigned  to  the  qk,nth  subcarrier,  where  qk,ne[  1,  ...,  M]. 
The  subcarrier  allocation  is  performed  by  a  control  unit  at  the  base  station  based  on  the  CSI 
received  from  the  mobiles,  and  the  channel  allocation  instructions  \qk,n}  are  sent  to  each  mobile 
via  a  forward  control  channel.  As  in  [5],  more  than  one  substream  of  the  same  user  can  hop  onto 
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the  same  subcarrier.  The  total  number  of  substreams  assigned  to  the  /nth  subcarrier  is  Pm.  Note 
the  bit  energy  Eb  is  the  same  for  each  substream  independent  of  its  subcarrier  allocation.  Thus, 
power  control  that  adapts  to  large-scale  fading  variation  is  assumed.  Better  system  performance 
can  be  achieved  by  combining  subcarrier  allocation  and  per-subcarrier  power  control,  i.e.,  the 
power  control  that  compensates  for  small-scale  fading  variations.  However,  the  latter  would 
increase  the  required  feedback  rate  from  the  base  station  to  the  mobile  users.  Moreover,  joint 
allocation  of  subcarrier  and  power  is  very  complex.  In  contrast,  our  method  exploits  both 
frequency  and  multiuser  diversity  while  maintaining  relatively  low  complexity. 

While  it  is  not  practical  to  maintain  orthogonality  among  substreams  of  different  users  [5,  11], 

T 

it  is  feasible  to  require  it  for  each  user  k,  i.e.,  cknCkM’=()  for  nVn.  These  orthogonal  codes  are 
obtained  by  first  assigning  random  spreading  codes  to  all  users  and  then  multiplying  the 
substreams  of  user  k  by  orthogonal  Walsh  codewords  [10,  18]. 

The  low-pass  equivalent  transmitted  signal  for  the  Mi  user  is  the  sum  of  modulated  substream 

N 

signals  Sk(t)=Y.  xk  „(t)cxp(-j2  n/'  t) ,  where/,  is  the  subcarrier  frequency  offset  from  the 

k,n 


carrier  frequency  f  for  the  qk,n th  subcarrier. 

B.  Channel  Model  and  Imperfect  CSI 

We  assume  slowly  varying  flat  Rayleigh  fading  for  each  subcarrier.  The  channel  coefficients 
Yk,m(t)  of  all  subcarriers  are  identically  distributed  with  E[\ykm(t)\  ]=1,  and  are  independent  for 
different  users.  Denote  the  amplitude  akim(t)=\/k,m(t)\  and  the  phase  (f>km(t)=Zykm(t). 

To  assess  realistic  system  performance,  we  employ  correlated  subcarriers.  Assuming  the 
propagation  delay  is  exponentially  distributed  with  rms  delay  spread  <Jd  [19,  p.  50],  the  cross¬ 
correlation  between  y kyin(i)  and  y /,.,„(/)  with  index  difference  A m=m-m'  is  given  by  [20] 


rfAmAf)  = 


E 

l 

5 

*  -ttf 
r*” 

i _ i 

l+2jnAmAf(Jd 

VE[lY*,m(0|2]  E[|y k,m 

-<0|2]  l+(27tA/?zA/arf)2 

(2) 
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In  the  reverse  link,  the  base  station  accurately  estimates  the  current  CSI  of  all  users  at  each 
subcarrier  by  employing  pilot  signals  [21].  However,  due  to  the  delay  caused  by  the  calculation 
of  allocation  variables  and  feeding  the  allocation  assignments  back  to  mobile  users,  the  CSI  used 
for  allocation  is  different  from  the  actual  channel  conditions  experienced  during  the 
transmission.  An  efficient  approach  to  reliable  transmission  in  the  presence  of  CSI  mismatch  is 
to  predict  future  CSI  Ykm(t)  based  on  outdated  observations  using  the  autoregressive  (AR)  model- 
based  linear  prediction  (LP)  [14,  22],  Since  the  random  variables  and  Ykm(t)  are 

jointly  Gaussian,  the  reliability  of  prediction  is  characterized  by  the  subcarrier- independent 

cross-correlation  rj  = E  [  %„,  ( t) A,  (0  W  E  [ Yk m  ( 0 1 2  ]  E  [  | Yk (0 1 2 ]  [23,24].  The  normalized  prediction 
mean  square  error  (NPMSE)  is  given  by  l-r\  [23],  For  example,  r|=0.925,  0.95,  0.975  and  1.0 
correspond  to  the  NPMSE=0.144,  0.098,  0.049,  and  0,  i.e.,  poor,  degraded,  reliable  and  perfect 
prediction,  respectively  [22,  23].  Implementation  of  reliable  practical  prediction  algorithm  for 
MC-CDMA  with  AFH  is  outside  the  scope  of  this  paper.  To  evaluate  performance  for  prediction 
accuracy  specified  by  the  reliability  parameter  r)  in  simulations,  we  generate  predicted  channel 
coefficients  as  Ykm(t)=  rj  Yk,m(t)+co(t),  where  ait)  is  white  Gaussian  noise  with  variance  ) 

that  is  independent  of  Ykm(t)- 
C.  The  Receiver 

At  the  base  station,  the  low-pass  equivalent  received  signal  is 

K  N 

K0=Z  Z  a.  fiK«(7-T/c)exp  SjX (phn{t)-2Kfqk  f\ }  +  ffi)  (3) 

k=\  n= 1  ’  Kn  K  J 

where  v{t)  is  the  complex  additive  white  Gaussian  noise  with  power  spectral  density  N0  ,  ik  is  the 
timing  offset  of  the  Alh  user,  and  (pu^(t)=2nfq kJTk+  <f>k,q kn(t)  ■  In  this  paper,  we  assume  perfect 

knowledge  of  the  timing  and  frequency  offsets  for  all  users,  and  perfect  CSI  at  the  receiver.  As 
discussed  in  [1 1],  we  adopt  the  quasi-synchronous  assumption  that  is  suitable  for  the  reverse  link 
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in  practice,  i.e.,  the  timing  of  all  users  is  aligned  within  a  small  synchronization  window  of  one 
to  a  few  chips  in  length. 

Without  loss  of  generality,  we  focus  on  detection  of  the  symbol  buM)  at  the  output  of  the 
correlator  at  the  receiver  and  drop  the  symbol  index  /.  Denote  /k,m=Tk,m(t) ,  ak,m=0Ck,m(t) ,  and 

<Pk,m=<Pk,m(t)  for  all  k  and  m  during  this  symbol  period. 

Define  U„,  of  size  Pm  as  the  index  set  of  all  substreams  allocated  to  the  /nth  subcarrier,  i.e. 
U m={(k,  n)  I  qk,n=m}.  Without  loss  of  generality,  we  replace  the  substream  index  pair  (k,  n)eUm 
with  a  single  index  p,  and  assume  the  corresponding  0,  t /,=().  and  ap=a^m.  The  output  of  the 

correlator  for  the  pi\\  substream  is  given  by  [5,  6,  8] 

Zp=Sp+Ip+Np,  pe\\,...,Pm\  (4) 

where  Sp=\jEi,apbp  carries  the  desired  bit  information,  and  Ip  is  the  interference  term 

~IP=^Eb  Z  apbp'pp<p  (5) 

p's  u,„ 

&p'*p 

where  pp'p  is  the  cross-correlation  between  the  waveforms  of  corresponding  substreams  [12] 

Tb  PG- 1  PG- 1 

P{k',n%k,n)=VLV\j((Pk,,m-(Pk,m)\  f  Z  [ck',n{h)h{t~i\Tc-Xk)\  Z  [ck,n(h)h(t-i2Tc)]dt  (6) 

q  /i=0  i’r= 0 

In  (5),  the  inter-user  interference  from  adjacent  symbols  and  other  subcarriers  is  negligible  [5,  8, 
10,33],  Thus,  the  interference  term  in  (5)  is  due  primarily  to  the  symbols  of  other  users 
transmitted  on  the  mth  subcarrier  that  start  in  the  same  synchronization  window  as  the  desired 

symbol  bp.  Finally,  the  covariance  between  the  zero-mean  Gaussian  noise  tenns  Np  in  (4)  is 
E(NpNP')=No  pP> 

D.  SINR  Analysis 

The  received  signal  (4)  is  detected  and  demultiplexed  at  the  receiver  to  make  decisions  (see 
Fig.  1).  The  analytical  SINR  of  the  detector  is  utilized  by  the  allocation  method.  The  SINR 
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k(m,p)  for  the  /?th  substream  of  the  mth  subcarrier  is  given  by  the  ratio  of  the  desired  signal 
power  to  the  total  power  of  residual  interference  and  noise  at  the  input  to  the  decision  device. 
The  SINR  of  the  matched  filter  (MF)  detector  is  [5,  12] 

E  1,0.1 

AMf,(m,p)  =  Pm  (7) 

~PqEj  E(p\p)a2p'+N0 

p- i 

where 

A(p',p)AA[(f,„'),(«:,»)]^;  ™ise  (8) 

In  (7),  we  assume  perfect  synchronization,  since  this  assumption  results  in  the  worst  case  system 
performance  [10,17,18,33]. 

When  the  linear  decorrelating  detector  (LDD)  is  utilized,  MAI  is  completely  eliminated  in  (4), 
so  the  SINR  is  given  by  the  SNR  [12,  p.  249] 

P  h  (J. 

ALDD{m,p)- N^)pp  (9) 

where  the  PmxPm  cross-correlation  matrix  R  has  components  Kp>p<=ppp<  (see  (6)). 

The  SINR  of  the  linear  minimum  mean-square  error  (MMSE)  detector  can  be  computed  for  the 
/?th  substream  as  [12,p.  293-298,  25] 

AMMSE{m ,/?)=j  (l+^v^“)  }  0°) 

where  A=diag{a!,  a2,  . .  .,a PJ,  and  I  is  kxk  identity  matrix. 

The  simple  successive  interference  canceller  (SIC)  is  also  employed  in  this  paper.  Without  loss 
of  generality,  assume  that  ai>a2>  ...>a ,Pm  [12],  The  SINR  of  the  SIC  can  be  approximated  as  [12, 


7 


where  Pe,  sic(m  ,p)~  Q  [V  2  A  sidjn ,/;)]  is  the  BER  estimate  for  the  /;th  substream, 

00 

Q(x)=fexp(-t2/2)dt/'\[2n,  and  A (p\p)  is  defined  in  (8). 

III.  ALLOCATION  ALGORITHMS 

In  this  section,  we  discuss  the  limitations  of  previously  proposed  WF  allocation  algorithm  and 
describe  a  novel  allocation  method  for  the  reverse  link.  As  discussed  in  section  II  (see  Fig.  1), 
these  algorithms  are  implemented  at  the  base  station,  and  the  resulting  allocation  assignments  are 
sent  back  to  the  mobile,  where  they  are  employed  to  select  the  reverse  link  subcarriers  for  each 
substream.  Thus,  the  computational  load  associated  with  channel  allocation  is  borne  entirely  by 
the  base  station,  which  can  afford  much  greater  complexity  than  the  mobiles.  Moreover,  the 
feedback  load  is  low  since  only  a  few  allocation  bits  need  to  be  fed  back  at  a  low  rate  on  the 
order  of  the  maximum  Doppler  shift. 

All  allocation  methods  employ  the  SINR  minimization  criterion  and  utilize  the  SINR 
expression  k(/?z,/?)  in  (7),  (9),  (10),  or  (1 1)  that  depends  on  the  type  of  the  receiver  utilized  at  the 
base  station.  Otherwise,  the  allocation  methods  do  not  depend  on  the  receiver  design  and  can  be 
easily  adapted  to  various  receiver  choices. 

Subcarrier  allocation  methods  were  also  addressed  in  [26,  27,  28,  29,  30].  However,  these 
papers  focus  on  the  orthogonal  frequency-division  multiple  access  (OFDMA)  system,  and  the 
throughput  and  fairness  among  users  are  primary  concerns,  since  the  target  application  is  the 
internet  protocol  (IP)-based  packet  transmission  system.  In  contrast,  we  consider  MC-CDMA 
with  AFH  that  has  fixed  data  rate  for  all  users  in  the  system,  which  is  suitable  for  such 
applications  as  voice  and  circuit-switched  data  transmission.  Thus,  it  is  reasonable  to  use  the 
minimum  SINR  as  performance  measure  to  guarantee  the  quality  of  service  for  all  users. 
Moreover,  since  the  multi-access  interference  (MAI)  is  not  significant  in  OFDMA,  design  of  the 
allocation  algorithm  is  simpler  for  this  system  than  for  MC-CDMA. 


A.  The  Waterfilling  (WF)  Algorithm 

The  WF  algorithm  was  proposed  in  [5]  for  the  MF  receiver.  Note  that  this  method  is  not  the 
water-filling  solution  that  achieves  the  information-theoretic  capacity,  although  its  design  was 
inspired  by  the  latter.  The  WF  algorithm  is  an  iterative  allocation  algorithm,  i.e.,  one  substream 
is  allocated  in  each  iteration.  The  sequence  of  KN  iterations  is  divided  into  N  consecutive  groups, 
with  K  iterations  per  group.  In  each  group,  one  substream  per  user  is  allocated.  During  each 
iteration  of  the  nth  group,  all  users  whose  substreams  have  not  yet  been  allocated  in  this  group 
select  the  subcarrier  with  the  highest  SINR,  and  the  user  that  has  the  lowest  selected  SINR  level 
is  assigned  to  that  subcarrier  [5].  The  protection  of  weak  users  offered  by  the  order  of  subcarrier 
allocation  in  the  WF  method  is  limited  due  to  two  reasons.  First,  weak  users  already  assigned  to 
given  subcarrier  cannot  contribute  large  interference  level  that  would  prevent  allocation  of  a 
much  stronger  user’s  substream  to  this  subcarrier.  Second,  the  WF  algorithm  does  not  take  into 
account  the  MAI  impact  of  assigning  a  new  substream  on  the  substreams  already  allocated  to 
certain  subcarrier.  As  a  result,  the  BER  of  the  WF  allocation  algorithm  is  very  high  for  some 
realizations  of  channel  fading  coefficients.  This  case  is  demonstrated  in  the  following  example. 


Example  Consider  a  system  with  K=N=M=2,  PG=  16,  Etr :1,  Ao=0.01  and  r|=1.0.  Suppose  fading 


2  2  2  2 

amplitudes  are  represented  by  {ai  i= 0.1,  cq  r=0.09,  a 2,i=l,  a2, 2=0.1}. 

The  WF  algorithm  [5]  performs  the  allocation  in  four  iterations.  The  resulting  allocation  is 
{gi,i=l,  q i,2=2,  <72,1=1,  <72,2=1 } .  Suppose  that  substream  (k=  1,  n=  1)  corresponds  to  substream p=\ 
of  the  first  subcarrier.  From  (7),  the  BER  of  this  substream  is  approximated  as 


Pe,Mfim=q  i ,  i 1 = 1 ,  P = 1  )*Ql 


I  2x0,1 
0.01+2/16 


*0.11. 


A  better  allocation  solution  is  to  assign  the  substreams  of  the  first  user  to  the  first  subcarrier 
and  the  substreams  of  the  second  user  to  the  second  subcarrier,  i.  e.,  {<71,1=1,  <71,2=1,  <72,1=2, 

<72,2=2}.  The  resulting  BERs  are  Pe>MA.tn,  p)~Q(fj  2x0. 1/0.0 1)«3. 8xl0'6  for  m,pe{  1,  2}. 
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This  example  illustrates  that  the  WF  algorithm  has  poor  performance  when  the  instantaneous 
fading  amplitude  of  one  user  (user  2  in  this  example)  is  much  larger  than  that  of  the  other  user. 
While  such  events  are  unlikely,  they  dominate  the  average  BER  in  moderate  to  high  SNR  region. 
Therefore,  it  is  desirable  to  design  an  allocation  algorithm  that  overcomes  the  limitations  of  the 
WF  method. 

B.  Proposed  Allocation  Algorithm 

Our  allocation  method  has  the  following  optimization  criterion.  As  in  [5],  we  focus  on  the 
performance  of  the  substream  with  the  lowest  SINR  since  the  error  events  associated  with  this 
substream  dominate  the  error  rate.  The  optimization  objective  is  to  maximize  the  SINR  of  the 
allocated  substream  with  the  lowest  SINR,  i.e.,  to  find  the  set  of  allocation  variables  that 

(90=arg  maxj  min  k(m,p)  \  (12) 

n  ]  '»e[l . M\  i 

~  l  pe[  1,. J 

where  Q={qk,n\ke[l,...,K\,  ne[l,...,N]}  is  the  set  of  all  possible  allocation  variables,  and  h(m,p) 
is  the  SINR  calculated  from  (7),  (9),  (10),  or  (11),  depending  on  the  type  of  the  detector 
employed  at  the  receiver.  In  (12),  SINRs  of  all  substreams  are  used  for  allocation,  so  the  weak 
user  is  protected  well.  Thus,  better  perfonnance  is  expected  relative  to  the  WF  algorithm. 

Theoretically,  the  search  over  the  elements  of  Q  produces  the  optimal  solution  for  (12).  There 
are  (A+M-1)!/[(M-1)!A!]a  possible  ways  to  allocate  N  equivalent  substreams/user  into  M 
distinct  subcarriers  [33].  Thus,  the  computational  complexity  of  finding  the  exact  solution  to  (12) 
is  very  high,  and  a  suboptimal  method  with  moderate  complexity  is  desirable. 

Inspired  by  the  iterative  algorithms  in  [5,9,26,27],  we  design  an  iterative  allocation  method  to 
find  a  suboptimal  solution  for  (12).  In  this  algorithm,  one  substream  is  assigned  at  each  iteration, 
and  all  substreams  are  allocated  after  KN  iterations.  For  simplicity,  the  substreams  are  assigned 
consecutively,  i.e.  the  substream  ( k ,  n)  is  allocated  at  the  [K(n  -l)+k]th  iteration.  It  is  shown  in 
[33]  that  this  simple  assignment  order  results  in  negligible  performance  degradation  relative  to 
more  sophisticated  methods.  The  algorithm  is  based  on  the  following  idea.  Suppose  the 
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substream  ( k ,  n)  is  allocated  to  the  subcarrier  m.  We  observe  that  in  this  case  only  the  substream 
( k ,  n)  and  the  substreams  previously  assigned  to  this  subcarrier  experience  the  SINR  degradation. 
Denote  the  set  of  such  substreams  U'm.  For  each  m  in  we  determine  the  lowest  SINR 

level  in  U'm  assuming  that  the  substream  (k,n)  is  allocated  to  subcarrier  m.  Then  find  the 
subcarrier  m0  that  produces  the  maximum  value  of  this  lowest  SINR  among  all  M  values  and 
allocate  the  substream  (k,n)  to  subcarrier  m0.  In  our  allocation  process,  the  substreams  are 
allocated  consecutively,  and  the  interference  from  a  subset  of  all  substreams  is  taken  into  account 
for  each  allocation.  Therefore,  it  results  in  a  suboptimal  solution  to  (12).  However,  it  improves 
upon  the  WF  method  in  [5]  since  it  protects  previously  allocated  weaker  users.  The  proposed 
improved  allocation  algorithm  consists  of  the  following  steps: 

SO)  Initialize  Pm= 0  and  Um={0}  for  all  m  in  [1 M\ 

for  n=  1  to  N 
for  A=1  to  K 

51)  Augment  U'm={Um,  (k,n)\  for  me[  1,  ...,  M ].  Let  k(m,p)  be  the  SINR  of  the  /;th 
substream  (pe[l,...,Pm+l])  assuming  that  substream  (k,n)  is  assigned  to  the  mth 
subcarrier. 

52)  Find  m0  that  satisfies 

/«0=arg  max  j  min  A (m,p)  \ 

53)  Assign  substream  (k,n)  to  the  m0th  subcarrier,  i.e.  set  qk,n=m0.  Then  update 
P  —P  +1  anHTT  =TT' 

1  m0  1  m0  1 5  ai  1U  um0  ^  m0* 

When  perfect  CSI  is  available  at  the  base  station,  the  SINR  in  step  SI  of  the  improved 
allocation  algorithm  is  calculated  using  the  expression  (7),  (9),  (10),  or  (11),  depending  on  the 
type  of  the  detector  employed  at  the  receiver.  As  discussed  earlier,  this  perfect  knowledge  of  CSI 
assumed  in  the  calculation  of  X(m,  p)  is  usually  unavailable.  Since  an  algorithm  that  takes  the 
prediction  error  into  account  is  very  complex,  we  focus  on  a  simpler  approach  that  employs  the 
predicted  fading  amplitude  ap  instead  of  ap  in  the  SINR  ),(m,p)  and  obtains  the  SINR  estimate 
%(m,p).  Thus,  the  predicted  CSI  is  treated  as  if  it  was  perfect,  and  Mm,p)  is  used  as  performance 
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measure  in  step  SI  of  our  allocation  method.  The  sensitivity  of  allocation  algorithm  to  the 
prediction  errors  is  tested  using  simulations  in  section  V. 

C.  Complexity  Reduction 

In  our  allocation  algorithm,  the  computational  complexity  is  mainly  detennined  by  the 
calculation  of  the  SINR.  For  MF  and  SIC  receivers,  the  calculation  of  the  SINR  given  by  (7)  and 
(11),  respectively,  involves  only  simple  scalar  operations.  On  the  other  hand,  the  allocation 
algorithms  for  the  LDD  or  MMSE  detector  are  complex  due  to  the  matrix  inversion  in  (9)  and 
(10),  which  has  to  be  performed  for  every  U'm  in  step  SI  of  the  improved  allocation  algorithm. 
However,  it  is  possible  to  reduce  complexity  by  recursive  computation  of  this  inverse  as 
discussed  below  for  the  case  of  the  LDD  at  the  receiver.  In  this  recursion,  for  each  m  in  step  SI, 
denote  the  PmxPm  cross-correlation  matrix  of  the  substreams  previously  assigned  to  subcarrier  in 
as  RfPm).  The  inverse  of  this  matrix  [RfP,„)]  1  has  been  computed  during  an  earlier  iteration.  In 
the  current  iteration,  we  need  to  evaluate  the  SNR  XLDD(m,p),  pe\  1,  ...,  P,„+ 1]  when  assignment 
of  new  (Rm+l)st  substream  is  examined  for  subcarrier  m  in  step  SI.  This  computation  requires 
the  inversion  of  an  augmented  correlation  matrix 


R(P,n+l)= 


R (Pm)  P 

P"  1 


(13) 


where  p=[pP  +1 ,,  pPm+1  2,. . .,  p/:  _ ,  r]r  is  the  vector  of  the  cross-correlation  coefficients  between 
the  new  substream  and  Pm  existing  substreams  determined  by  (6).  The  inverse  of  (13)  can  be 
computed  by  exploiting  the  existing  [R(/f„)]  1  as  follows  [31,  eq.  (35)]: 


[R(  Pm+\  )]  = 


[R  {Pm)f+  p  [R(Pm)YlppH[R(Pm)Tl  -p[R(Rm)r'p 


L-pPH[R(Rm)r 


p 


(14) 


where  p=[l-p//[R(/,„;)]  'p]  This  recursion  requires  about  1. 5 (Pm+1)2  complex  multiplications 
and  divisions  [33],  which  is  significantly  less  than  the  (Pm+  l)log27  operations  of  efficient 
Strassen’s  method  for  direct  matrix  inversion  [32].  The  SINR  calculation  for  MMSE  detector 
(10)  can  be  simplified  similarly  [33]. 
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IV.  NUMERICAL  RESULTS 


First,  we  compare  the  average  BER  (over  all  users  and  fading  realizations)  of  the  WF 
algorithm,  the  optimal  allocation  method  (12)  and  the  proposed  algorithm  for  the  MF  detector  in 
Fig.  2.  As  in  [5],  the  assumptions  of  perfect  CSI  knowledge  and  independent  fading  subcarriers 
are  adopted  in  this  simulation,  as  well  as  in  Fig.  3  and  4.  Due  to  very  high  complexity  of  the 
optimal  allocation,  small  number  of  substreams  is  employed  in  this  comparison.  We  observe  that 
the  proposed  algorithm  incurs  a  small  penalty  relative  to  the  optimal  allocation  method  for  high 
SNR  and  outperfonns  the  latter  for  low  SNR  since  the  average  BER  is  not  the  optimality 
criterion.  As  expected,  the  WF  algorithm  has  the  worst  perfonnance  among  the  three  allocation 
algorithms  in  moderate  to  high  SNR.  Fig.  2  also  demonstrates  that  the  AFH  system  with  the 
proposed  allocation  algorithm  provides  2-3dB  power  gain  over  the  non-adaptive  MC  DS-CDMA 
with  MRC  [8]  in  this  multiuser  scenario. 

Fig.  3  compares  the  average  BER  of  the  WF  method  and  the  proposed  allocation  algorithm  for 
the  MF  receiver  and  a  larger  number  of  users  in  the  system.  We  observe  that  the  perfonnance 
gain  of  the  proposed  algorithm  over  the  WF  method  increases  as  the  SNR  grows.  We  have  also 
investigated  the  modification  of  the  WF  algorithm  when  the  LDD  is  employed  in  the  receiver  in 
[33].  Significant  loss  was  also  observed  relative  to  the  proposed  allocation  method  in  this  case. 
The  comparison  in  Fig.  3  corroborates  our  discussion  on  the  limitations  of  the  WF  algorithm  in 
section  III-A. 

Fig.  4  shows  the  average  BER  of  the  proposed  allocation  algorithm  for  moderate  system  load 
(. K  comparable  to  M)  combined  with  receivers  discussed  in  section  II.  To  illustrate  the  potential 
of  combined  subcarrier  allocation  and  MUD,  the  single  user  bound  (SUB),  i.e.,  the  BER  of  the 
system  with  only  one  user,  is  included  in  the  figure.  The  linear  methods  (the  MMSE  detector  and 
LDD)  have  much  better  performance  than  the  SIC  and  the  MF  detectors  for  moderate  to  high 
SNR.  The  linear  detectors  approach  the  SUB  in  this  system  because  the  adaptive  allocation 
algorithm  separates  the  substreams  with  high  cross-correlations  into  different  subcarriers.  This 
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adaptive  allocation  policy  also  reduces  the  effect  of  noise  enhancement  that  degrades  LDD  in 
conventional  CDMA  channels,  so  it  is  only  slightly  outperformed  by  the  MMSE  detector  in  our 
system.  We  also  note  that  when  the  allocation  method  is  employed  with  the  SIC,  it  is  impaired  in 
the  first  few  stages  by  the  inaccuracy  of  the  SINR  approximation  (11)  when  the  number  of 
substreams  is  small  [12]. 

We  investigate  the  effect  of  imperfect  CSI  in  Fig.  5.  Observe  that  while  the  BER  of  the 
adaptive  system  with  LDD  increases  as  the  reliability  of  CSI  degrades,  it  is  still  smaller  than  the 
SUB  of  non-adaptive  MC  DS-CDMA  even  for  r|=0.95,  which  is  easily  achievable  in  realistic 
scenarios  [15,  22]. 

Fig.  6  illustrates  the  impact  of  correlated  subcarriers  on  the  BER  of  the  adaptive  system  when 
perfect  knowledge  of  the  CSI  is  assumed.  As  expected,  the  performance  degrades  as  A faj 
decreases.  However,  for  A fad  as  small  as  0.3,  the  performance  is  very  close  to  the  case  of 
independent  fading  subcarriers  (A/ay=oo).  Observe  that  the  BER  gap  between  the  LDD  and  SIC 
narrows  as  A fod  decreases  in  Fig.  6  since  fading,  not  MAI,  becomes  the  dominant  source  of 
performance  degradation  as  diversity  reduces  for  small  A foc/.  Similarly,  it  is  shown  in  [33]  that 
performance  gap  between  these  two  detectors  narrows  when  the  predicted  CSI  is  unreliable  in 
the  adaptive  system.  The  complexity  comparison  in  [33]  reveals  that  the  allocation  algorithms  for 
SIC  and  linear  detectors  have  similar  complexity.  Therefore,  the  LDD  or  MMSE  detectors 
provide  the  best  complexity/performance/robustness  trade-off  in  this  adaptive  system. 

In  Fig.  7  and  8,  we  investigate  a  practical  system  design  with  the  available  bandwidth  set  to 
1.25MHz,  the  bandwidth  of  single  carrier  in  IS-95  and  CDMA2000  standard.  The  nns  delay 
spread  of  the  channel  is  assumed  acf=  1  us,  which  is  typical  for  the  urban  environment  [18,  p.  200]. 
Since  the  coherence  bandwidth  is  approximately  1/5gj=200KHz  [18,  p.  202],  we  divide  the 
available  bandwidth  into  12  subcarriers  with  the  subcarrier  bandwidth  of  102.4KHz  to  assure  flat 
fading  for  each  subcarrier.  Assuming  N=  8  substreams  and  PG=6 4,  the  system  supports  a 
12.8Kbps  data  service  to  each  user.  Moreover,  for  comparison,  we  have  simulated  a  non- 
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adaptive  MC  DS-CDMA  system  with  the  same  system  bandwidth,  number  of  subcarriers,  and 
data  rate  as  for  the  adaptive  system.  To  maintain  the  data  rate,  the  time-domain  processing  gain 
of  eight  is  required  in  the  latter  system.  While  in  the  original  MC  DS-CDMA  design,  the  same 
spreading  code  was  used  for  all  subcarriers  of  one  user  [8],  we  employ  different  random 
spreading  codes  for  different  subcarriers  as  proposed  in  [13]  to  reduce  the  variance  of  the  cross¬ 
correlation  between  the  waveforms  of  different  users. 

Fig.  7  shows  the  BER  of  the  adaptive  system  for  varying  numbers  of  users.  For  all  detectors, 
the  BER  level  is  maintained  provided  the  number  of  users  K  is  sufficiently  small.  This  is  due  to 
the  fact  that  the  fading  induced  near-far  problem  can  be  easily  suppressed  for  a  small  number  of 
users.  While  the  BER  of  the  MF  receiver  degrades  for  K>M,  the  MUDs  maintain  reliable 
performance  unless  the  system  load  is  very  large,  i.  e.,  when  K»M.  The  MMSE  detector  has  the 
best  performance,  but  its  gain  over  the  LDD  is  very  small  as  discussed  above. 

In  Fig.  8,  the  BERs  of  the  adaptive  and  non-adaptive  systems  with  the  LDD  are  compared.  For 
the  non-adaptive  system,  this  MUD  is  derived  from  the  RAKE  LDD  designed  for  the  multipath 
fading  channel  in  [34],  Even  with  poorly  predicted  CSI  (r|=0.925),  the  AFH  system  outperforms 
the  non-adaptive  system.  For  target  BER  of  10'  ,  the  adaptive  system  can  support  from  12  to  63 
users,  depending  on  the  quality  of  the  predicted  CSI,  while  only  two  users  are  feasible  for  the 
non-adaptive  system. 

In  previous  simulations,  the  power  control  that  adapts  perfectly  to  large-scale  fading  is 
assumed,  i.e.,  Eb  is  the  same  for  all  users.  To  investigate  the  effect  of  the  near-far  problem  caused 
by  imperfect  power  control  on  the  proposed  allocation  algorithm,  the  SNR  of  randomly  chosen 
desired  user  is  set  to  lOdB,  while  K-l  interfering  users  have  SNR=20dB.  Given  the  same  system 
parameters  as  in  Fig.  7,  the  simulated  BER  of  the  desired  user  is  shown  in  Fig.  9.  Note  that 
MUDs  in  this  figure  have  lower  BER  than  for  the  perfect  power  control  case  in  Fig.  7,  and  the 
BERs  of  linear  detectors  are  very  close  to  the  system  SUB  even  for  heavily  loaded  near-far 
scenario.  The  reason  is  that  the  weakest  user’s  substreams  are  assigned  into  their  favorite 


15 


subcarriers  by  our  allocation  algorithm,  while  the  interfering  stronger  users  are  competing  for 
remaining  subcarriers.  In  previous  figures,  our  allocation  policy  was  shown  to  alleviate  small- 
scale  fading-induced  near-far  problem.  Fig.  9  demonstrates  that  combined  proposed  allocation 
algorithm  and  MUD  are  also  very  effective  in  compensating  for  large-scale  fading  and  imperfect 
power  control. 

V.  CONCLUSIONS 

A  novel  allocation  algorithm  that  outperforms  the  WF  algorithm  for  MC-CDMA  system  with 
AFH  for  the  reverse  link  was  proposed.  The  proposed  allocation  algorithm  is  combined  with 
MUDs  using  the  SINR  analysis.  It  is  demonstrated  that  combined  improved  subcarrier  allocation 
and  MUD  is  very  effective  in  mitigating  MAI,  resulting  in  much  larger  system  capacity  than  for 
the  non-adaptive  MC  DS-CDMA  system.  This  conclusion  is  confirmed  for  realistic  rapidly 
varying  fading  channels  with  correlated  subcarriers,  CSI  mismatch,  and  imperfect  power  control. 
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Fig.  1  System  diagram  of  MC-CDMA  with  AFH  for  reverse  link 
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Fig.  2  BER  performance  of  adaptive  MC-CDMA  with  WF,  optimal  and  proposed  allocation 
algorithms  (N=l  substream/user,  PG=64  and  perfect  CSI  knowledge),  and  non-adaptive  MC  DS- 
CDMA  with  MRC  (time  domain  processing  gain  PG=64).  For  both  systems,  K=  6  users,  M=  4 
subcarriers,  independent  fading  subcarriers  (A/g</=oo),  and  MF  receiver  are  assumed. 


Fig.  3  BER  performance  of  adaptive  MC-CDMA  with  MF  receiver  for  K=  16,  M=  8,  PG= 64, 

Afod=co,  and  r)=l  .0 
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Fig.  4  BER  performance  of  MC-CDMA  with  AFH  for  K=\6,  N=  8,  and  r|=1.0.  In  all  systems, 

M=  8,  PG= 64,  and  A/a(/=cc. 


Fig.  5  Comparison  of  SUB  of  non-adaptive  MC  DS-CDMA  with  the  BER  of  MC-CDMA  with 
AFH  and  LDD  for  K=\6,  N=  8,  and  varying  CSI  reliability  p.  In  all  systems,  PG=64,  M=  8,  and 

Afoj=cc. 


21 


Fig.  6  BER  performance  of  MC-CDMA  with  AFH  and  LDD/SIC  for  K=  16,  M=8,  N=8,  PG= 64,  and  perfect 

CSI  (ri=1 .0)  at  the  allocation  algorithm. 


Fig.  7  BER  vs.  number  of  users  in  the  adaptive  system  with  improved  allocation  algorithm  for 
SNR=10dB,  M=12,  N=S,  PG= 64,  r|=0.95,  and  AfarO.1024 
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Fig.  8  BER  vs.  number  of  users  in  the  system  with  LDD  for  SNR=10dB,  M=  12,  7V=8,  PG= 64, 

and  A/arf=0.1024. 


Fig.  9  BER  vs.  number  of  users  in  the  adaptive  system  with  imperfect  power  control  for  /W=12,  A/=8, 
PG=64,  ri=0.95,  and  Afad-0. 1 024.  The  SNR  is  1 0dB  for  desired  user  and  20dB  for  other  users. 
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Table  I  Notation 


K 

number  of  users 

M 

number  of  subcarriers 

N 

number  of  substreams/user 

k,  k' 

user  indices 

m,  m ' 

subcarrier  indices 

n,  n ' 

substream  indices  for  a  given  user 

P>  P' 

substreams  indices  for  a  given  subcarrier 

( k,n ) 

the  nth  substream  of  the  kth  user 

ik 

timing  offset  of  the  Mi  user 

Eb 

bit  energy 

No 

power  spectral  density  of  channel  Gaussian  noise 

Tc 

chip  duration 

Tb 

bit  duration 

PG 

processing  gain=  Ti/Tc 

A  f 

subcarrier  bandwidth 

CT  d 

rms  delay  spread 

Qk,n 

allocation  variable  for  substream  (k,n) 

yk.m(t) 

fading  coefficient  on  the  mth  subcarrier  of  user  k 

A 

prediction  of 

A 

^k,mif)t  OL k,m(j ) 

A 

the  amplitudes  of  y and  Yk,m{t) 

$k,m(f) 

phase  of  y k.m(t) 

11 

A 

cross-correlation  between  yk,m(t)  and  yk,m(t) 

U  m 

set  of  all  substreams  assigned  to  subcarrier  m 

P 

1  m 

number  of  substreams  in  subcarrier  m 

Mm,p) 

SINR  for  the  /2th  substream  on  subcarrier  m 
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